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1  Introduction 

The  core  subject  of  my  research  supported  by  AFOSR  is  imaging  in  random  media  with  ap¬ 
plications  in  ultrasound  imaging  as  well  as  radar  remote  sensing  and  the  use  of  time-reversal  in 
array-to-array  communications. 

We  have  made  considerable  progress  in  both  imaging  and  time  reversal  during  the  three  year 
period.  In  addition  to  Borcea,  the  other  researchers  attached  to  this  project  were:  L.  Borcea,  J. 
Berryman,  C.  Tsogka,  A.  Kim  and  P.  Blomgren,  D.  Berebichez,  G.  Derveaux. 

Here  are  our  main  results 

•  Introduction  and  analysis  of  interferometric  array  imaging  for  cluttered  media.  Papers 

1,2, 3, 9. 

•  Analysis  of  Statisitcal  Stability  in  Time  Reversal.  Papers  4,10,11,19. 

•  Broadband  array  imaging  of  small  scatterers  in  clutter  using  the  singular  value  decomposition. 
Papers  5,6. 

•  Study  of  super-resolution  and  statistical  stability  in  Time  Reversal,  including  numerical  sim¬ 
ulations.  Papers  7,8. 

•  Time  reversal  communication  systems.  Analysis  of  cross-channel  interference  (spatial  focus¬ 
ing)  and  intersymbol  interference  (time  focusing).  Role  of  time-reversal  in  assessing  these 
issues.  Papers  13-17. 

•  Various  papers  (12,  18,  20)  in  wave  propagation  and  imaging  related  to  this  work  mathemat¬ 
ically. 
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IMAGING  and  TIME-REVERSAL 

1.  Interferometric  array  imaging  in  clutter,  Liliana  Borcea,  George  Papanicolaou  and 
Chrysoula  Tsogka.  Submitted  to  Inverse  Problems. 

Abstract: 

We  introduce  a  space-time  interferometric  array  imaging  functional  that  provides  statistically 
stable  images  in  cluttered  environments.  We  also  present  a  resolution  theory  for  this  imaging  func¬ 
tional  that  relates  the  space-time  coherence  of  the  data  to  the  range  and  cross-range  resolution  of 
the  image.  Extensive  numerical  simulations  illustrate  the  theory  and  address  some  implementation 
issues. 

2.  Theory  and  applications  of  time  reversal  and  interferometric  imaging,  Liliana 
Borcea,  George  Papanicolaou  and  Chrysoula  Tsogka.  Inverse  Problems,  vol  19,  (2003),  pp.  5139- 
5164. 

Abstract: 

In  time  reversal,  an  array  of  transducers  receives  the  signal  emitted  by  a  localized  source,  time 
reverses  it  and  re-emits  it  into  the  medium.  The  emitted  waves  back-propagate  to  the  source  and 
tend  to  focus  near  it.  In  a  homogeneous  medium,  the  cross-range  resolution  of  the  refocused  field 
at  the  source  location  is  Ao L/a,  where  Ao  is  the  carrier  wavelength,  L  is  the  range  and  a  is  the 
array  aperture.  In  a  noisy  (random)  medium,  the  cross-range  resolution  is  improved  beyond  the 
homogeneous  diffraction  limit  because  the  array  can  capture  waves  that  move  away  from  it  at  the 
source,  but  get  scattered  onto  it  by  the  inhomogeneities.  We  refer  to  this  phenomenon  as  super¬ 
resolution  of  the  time  reversal  process  in  random  media.  Super-resolution  implies  in  particular 
that,  because  of  multipathing,  the  array  appears  to  have  an  effective  aperture  that  is  greater 
than  a.  We  distinguish  two  types  of  effective  apertures:  (i)  the  narrow-band  one  ae  that  depends 
on  the  scattering  medium  and  (ii)  the  broad-band  one  Ae  which  is  determined  by  ae  and  by  the 
bandwidth  of  the  signal.  In  remote  sensing  regimes  (L  >  a),  the  influence  of  the  aperture  a  of  the 
array  on  the  focusing  resolution  is  negligible  and  focusing  is  essentially  determined  by  ae.  However, 
at  a  moderate  source  range  the  array  aperture  does  play  a  role  and  focusing  depends  on  Ae.  Since 
ae  depends  on  the  scattering  medium,  it  is  not  known.  In  this  paper  we  present  a  robust  procedure 
for  estimating  ae  from  the  signals  received  at  the  array.  The  interest  of  this  estimation  is  two¬ 
fold:  First,  knowing  ae  permits  assessing  quantitatively  super-resolution  in  time  reversal.  Thus,  it 
provides  an  estimation  of  the  refocusing  resolution  that  can  be  achieved  by  the  time  reversal  process 
in  cluttered  media  for  applications  such  as  spatially  localized,  secure  communications.  Second,  as 
we  show  in  this  paper,  ae  quantifies  in  an  explicit  way  the  loss  of  resolution  in  array  imaging. 

3.  A  resolution  study  for  imaging  and  time  reversal  in  random  media,  Liliana 
Borcea,  George  Papanicolaou  and  Chrysoula  Tsogka.  Contemporary  Mathematics,  Volume  333, 
G.  Alessandrini  and  G.  Uhlmann  editors,  AMS,  2003,  pp.  63-77. 

4.  Statistical  stability  in  time  reversal,  George  Papanicolaou,  Leonid  Ryzhik  and  Knut 
Solna.  SIAM  J.  on  Appl.  Math.,  64  (2004),  pp.  1133-1155. 

Abstract: 

When  a  signal  is  emitted  from  a  source,  recorded  by  an  array  of  transducers,  time  reversed  and 
re-emitted  into  the  medium,  it  will  refocus  approximately  on  the  source  location.  We  analyze  the 
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refocusing  resolution  in  a  high  frequency,  remote  sensing  regime,  and  show  that,  because  of  multiple 
scattering,  in  an  inhomogeneous  or  random  medium  it  can  improve  beyond  the  diffraction  limit. 
We  also  show  that  the  back-propagated  signal  from  a  spatially  localized  narrow-band  source  is 
self-averaging,  or  statistically  stable,  and  relate  this  to  the  self-averaging  properties  of  functionals 
of  the  Wigner  distribution  in  phase  space.  Time  reversal  from  spatially  distributed  sources  is 
self- averaging  only  for  broad-band  signals.  The  array  of  transducers  operates  in  a  remote-sensing 
regime  so  we  analyze  time  reversal  with  the  parabolic  or  paraxial  wave  equation. 

5.  Imaging  and  time  reversal  in  random  media,  Liliana  Borcea,  Chrysoula  Tsogka,  G. 
Papanicolaou  and  James  Berryman.  Inverse  Problems,  18  (2002),  pp.  1247-1279. 

Abstract: 

We  present  a  general  method  for  estimating  the  location  of  small,  well-separated  scatterers  in 
a  randomly  inhomogeneous  environment  using  an  active  sensor  array.  The  main  features  of  this 
method  are  (i)  an  arrival  time  analysis  of  the  echo  received  from  the  scatterers,  (ii)  a  singular  value 
decomposition  of  the  array  response  matrix  in  the  frequency  domain,  and  (iii)  the  construction  of 
an  objective  function  in  the  time  domain  that  is  statistically  stable  and  peaks  on  the  scatterers. 
By  statistically  stable  we  mean  here  that  the  objective  function  is  self-averaging  over  individual 
realizations  of  the  medium.  This  is  a  new  approach  to  array  imaging  that  is  motivated  by  time 
reversal  in  random  media,  analyzed  in  detail  previously.  It  combines  features  from  seismic  imaging 
like  arrival  time  analysis  with  frequency-domain  signal  subspace  methodology  like  Multiple  Signal 
Classification  (MUSIC).  We  illustrate  the  theory  with  numerical  simulations  for  ultrasound. 

6.  Statistically  stable  ultrasonic  imaging  in  random  media,  James  Berryman,  Liliana 
Borcea,  George  Papanicolaou  and  Chrysoula  Tsogka.  Journal  of  the  Acoustical  Society  of  America, 
112  (2002),  pp.  1509-1522. 

Abstract: 

Analysis  of  array  data  from  acoustic  scattering  in  a  random  medium  with  a  small  number  of 
isolated  targets  is  performed  in  order  to  image  and,  thereby,  localize  the  spatial  position  of  each 
target.  Because  the  host  medium  has  random  fluctuations  in  wave  speed,  the  background  medium 
is  itself  a  source  of  scattered  energy.  We  assume,  however,  that  the  targets  are  sufficiently  larger 
and/or  more  reflective  than  the  background  fluctuations  so  that  a  clear  distinction  can  be  made 
between  targets  and  background  scatterers.  In  our  numerical  simulations  we  use  nonreflective 
boundary  conditions  so  as  to  isolate  the  effects  of  the  host  randomness  from  those  of  the  spatial 
boundaries,  which  can  then  be  treated  in  a  separate  analysis.  We  show  that  the  key  to  successful 
imaging  is  finding  statistically  stable  functionals  of  the  data  whose  extreme  values  provide  estimates 
of  scatterer  locations.  The  best  ones  are  related  to  the  eigenfunctions  and  eigenvalues  of  the 
array  response  matrix,  just  as  one  might  expect  from  prior  work  on  array  data  processing  in 
complex  scattering  media  having  homogeneous  backgrounds.  The  specific  imaging  functionals 
studied  include  matched-field  processing  and  linear  subspace  methods,  such  as  MUSIC  (Multiple 
Signal  Classification).  But  statistical  stability  is  not  characteristic  of  the  frequency  domain,  which 
is  often  the  province  of  these  methods.  By  transforming  back  into  the  time  domain  after  first 
diagonalizing  the  array  data  in  the  frequency  domain,  we  can  take  advantage  of  both  the  time- 
domain  stability  and  the  frequency-domain  orthogonality  of  the  relevant  eigenfunctions. 

7.  Super-Resolution  in  Time-Reversal  Acoustics,  P.  Blomgren,  G.  Papanicolaou  and  H. 
Zhao.  Journal  of  the  Acoustical  Society  of  America,  Vol  111,  (2002),  pp.  230-248. 

Abstract: 

We  analyze  theoretically  and  with  numerical  simulations  the  phenomenon  of  super-resolution  in 
time-reversal  acoustics.  A  signal  that  is  recorded  and  then  re-transmitted  by  an  array  of  transduc¬ 
ers,  propagates  back  though  the  medium  and  refocuses  approximately  on  the  source  that  emitted 
it.  In  a  homogeneous  medium,  the  refocusing  resolution  of  the  time-reversed  signal  is  limited  by 


3 


* 


diffraction.  When  the  medium  has  random  inhomogeneities  the  resolution  of  the  refocused  signal 
can  in  some  circumstances  beat  the  diffraction  limit.  This  is  super-resolution.  We  give  a  theoretical 
treatment  of  this  phenomenon  and  present  numerical  simulations  which  confirm  the  theory. 

8.  Time  reversal  through  a  solid-liquid  interface  and  super-resolution,  Chrysoula 
Tsogka  and  G.  Papanicolaou.  Inverse  Problems,  18  (2002),  pp.  1639-1657. 

Abstract: 

We  present  numerical  computations  that  reproduce  the  time  reversal  experiments  of  Draeger, 
Cassereau  and  Fink,  where  ultrasound  elastic  waves  are  time  reversed  back  to  their  source  with 
a  Time  Reversal  Mirror  (TRM)  in  a  fluid  adjacent  to  the  solid.  We  also  show  numerically  that 
multipathing  caused  by  random  inhomogeneities  improves  the  focusing  of  the  back-propagated 
elastic  waves  beyond  the  diffraction  limit  seen  previously  in  acoustic  wave  propagation,  which  is 
called  super-resolution.  A  theoretical  explanation  of  the  robustness  of  super-resolution  along  with 
several  numerical  computations  that  support  this  explanation  is  in  our  previous  studies.  Time 
reversal  with  super-resolution  can  be  used  in  non-destructive  testing  and,  in  a  different  way,  in 
imaging  with  active  arrays. 

9.  Resolution  estimation  for  imaging  and  time  reversal  in  scattering  media,  Liliana 
Borcea,  G.  Papanicolaou  and  Chrysoula  Tsogka.  In  ” Mathematical  and  Numerical  Aspects  of 
Wave  Propagation  WAVES  2003”,  G.  Cohen,  E.  Heikkola,  P.  Joly  and  P.  Neittaanmaeki  Editors, 
Springer,  2003,  pp.  631-636. 

10.  The  parabolic  wave  approximation  and  time  reversal,  George  Papanicolaou,  Leonid 
Ryzhik  and  Knut  Solna.  Matematica  Contemporanea,  23,  (2002),  pp.  139-160. 

11.  Self-averaging  in  time  reversal  for  the  parabolic  wave  equation,  Guillaume  Bal, 
George  Papanicolaou  and  Leonid  Ryzhik.  Stochastics  and  Dynamics,  2,  (2002),  pp.  507-531. 

12.  Sensitivity  analysis  of  a  nonlinear  inversion  method  for  3D  electromagnetic 
imaging  in  anisotropic  media,  O.  Dorn,  H.  Bertete- Aguirre,  J.  Berryman  and  G.  Papanicolaou. 
Inverse  Problems,  18,  (2002),  pp.  285-317. 

TIME-REVERSAL  COMMUNICATIONS 

13.  Low  probabilty  of  intercept  and  intersymbol  interference  in  multiple-input 
single-output  time  reversal  communication  systems,  A.  Kim,  P.  Kyritsi,  P.  Blomgren  and 
G.  Papanicolaou.  Submitted  for  publication  2004. 

Abstract: 

We  study  a  multiple-input/single-output  underwater  communication  system  that  applies  time 
reversal  to  transmit  signals  so  that  they  focus  spatially  and  compress  temporally  on  the  intended 
receiver.  Our  simulations  model  an  underwater  acoustic  channel  as  a  waveguide.  The  results  verify 
the  theory  that  spatial  focusing  depends  strongly  on  the  product  of  the  delay  spread  and  the  band¬ 
width.  Because  of  spatial  focusing  at  the  intended  receiver,  this  system  has  a  low  probability  of 
intercept.  Even  though  signals  are  compressed  temporally,  there  remains  intersymbol  interference, 
especially  when  the  delay  spread-bandwidth  product  is  not  too  large.  However,  it  is  possible  to 
remove  nearly  all  intersymbol  interference  by  pre-equalizing  the  channel.  We  show  that  the  intro¬ 
duction  of  zero-forcing  pre-equalization  does  not  alter  significantly  the  spatial  focusing  properties 
of  time-reversal. 

14.  Predicted  Time  Reversal  Performance  in  Wireless  Communications  Using 
Channel  Measurements,  S.M.  Emami,  J.  Hansen,  A.D.  Kim,  G.  Papanicolaou,  A.J.  Paulraj,  D. 
Cheung  and  C.  Prettie.  IEEE  Comm.  Let.  2004. 

Abstract: 

Using  broadband  radio  wireless  measurements  in  an  indoor  environment  we  demonstrate  the 
remarkable  space-time  focusing  properties  of  signal  transmission  with  time  reversal. 

15.  Characterization  of  space-time  focusing  in  time  reversed  random  fields,  C. 
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Oestges,  A.D.  Kim,  G.  Papanicolaou,  A.J.  Paulraj.  IEEE  Trans.  Antennas  and  Prop,  vol.  53, 
(2005)  283-293. 

Abstract: 

This  paper  proposes  various  metrics  to  characterize  space-time  localization  resulting  from  ap¬ 
plication  of  time  reversal  techniques  in  richly  scattering  media.  The  concept  and  goals  of  time 
reversal  are  presented.  Then,  pertinent  metrics  describing  both  the  time  and  space  focalization 
effects  are  outlined.  Finally,  two  examples,  based  on  a  model  of  discrete  and  continuous  scattering 
media,  are  used  to  illustrate  how  the  proposed  metrics  vary  as  a  function  of  various  system  and 
channel  parameters,  such  as  the  bandwidth,  delay  and  angle  spreads,  number  of  antennas,  etc. 

16.  MISO  time  reversal  and  delay  spread  compression  for  FWA  channels  at  5GHz, 
P.  Kyritsi,  G.  Papanicolaou,  P.  Eggers  and  A.  Opera.  Antennas  and  Wireless  Propagation  Letters, 
3  (2004),  96-99. 

Abstract: 

Fixed  wireless  access  channels  in  the  5GHz  band  have  been  measured  with  8  element  uniform 
linear  antenna  arrays  at  both  ends,  thus  providing  an  8x8  configuration.  The  measurements  were 
performed  at  3  different  locations  in  downtown  Toronto.  The  application  of  the  time  reversal  (TR) 
technique  in  a  multiple  input-single  output  (MISO)  can  reduce  the  delay  spread  of  the  channel 
impulse  response  by  a  factor  of  2-3,  depending  on  the  power  allocation  scheme. 

17.  Application  of  Time-Reversal  with  MMSE  Equalizer  to  UWB  Communications, 
T.  Strohmer,  M.  Emami,  J.  Hansen,  G.  Papanicolaou  and  A.  Paulraj.  Global  Telecommunications 
Conference,  2004.  GLOBECOM  ’04.  IEEE  ,Vol  5,  29  Nov.-3  Dec.,  2004  Pages:3123  -  3127. 

OTHER  RELATED  WORK 

18.  Efficient  numerical  simulation  for  long  range  propagation,  K.  Huang,  G.  Papani¬ 
colaou,  K.  Solna,  C.  Thogka  and  H.  Zhao.  Submitted  for  publication. 

19.  Radiative  transport  limit  for  the  random  Schroedinger  equation,  G.  Bal,  G. 
Papanicolaou  and  L.  Ryzhik.  Nonlinearity,  vol  15,  (2002),  pp.  513-529. 

20.  High  frequency  behavior  of  the  focusing  nonlinear  Schroedinger  equation  with 
random  inhomogeneities,  A.  Fannjiang,  Shi  Jin  and  G.  Papanicolaou.  SIAM  Journal  on  Applied 
Mathematics,  vol  63,  (2003),  pp.  1328-1358. 
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